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A SIMPLE METHOD OF CALCULATING POWER-LAW VELOCITY 
PROFILE EXPONENTS FROM EXPERIMENTAL DATA 


By Jerry M. Allen 
Langley Research Center 


SUMMARY 

Analytical expressions for the effects of compressibility and heat 
transfer on laminar and turbulent shape factors H have been developed. 
Solving the turbulent equation for the power law velocity profile exponent 
N has resulted in a simple technique by which the N values of experimental 
turbulent profiles can be calculated directly from the integral 
parameters. Thus the data plotting, curve fitting, and slope measuring, 
which is the normal technique of obtaining experimental N values, is 
eliminated. The N values obtained by this method should be within the 
accuracy with which they could be measured. 


INTRODUCTION 


It has long been known that turbulent velocity profiles may be 
represented by a power law of the form 

(l) ^ ,(’> 

where N generally varies^from about 5 to 11 for turbulent velocity 
profiles. Figure l(a) shows one-fifth and one-eleventh power law profiles 

plotted in the conventional manner. 

Since the exponent N represents the general shape of the profile, 
the experimenter often obtains N values of his profiles as a useful 
index of profile shape. The usual method of obtaining these experimental 
N values is to plot the velocity profile in logarithmic form, fit 
a straight line through the data, and measure the slope of this line. 
Figure l{b) shows how these one-fifth and one-eleventh power law profiles 
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appear in logarithmic form. This process is rather tedious 
and time consuming, especially when a large number of profiles 

is involved. 

It is a common practice in the conduct of boundary layer experiments 
to obtain from the measured profiles certain integral parameters, such as 
displacement thickness 6* and momentum thickness 0 , which are descriptive 
of the character of the flow. These parameters are obtained by a process 
which is independent of the method normally used to define the profile N 
values. 

The objectives of this paper are to develop an analytical expression 
for the effects of compressibil ity and heat transfer on the turbulent shape 
factor, to compare this expression with a similar one for the laminar shape 
factor, and to demonstrate how this turbulent expression may be used as 
a simple method of calculating the N values of experimental turbulent 
velocity profiles, thus eliminating the need of logarithmic data plotting, 
curve fitting, and slope measuring. 

nomenclature 

H Shape factor, 6*/0 

M Mach number 

N power law velocity profile exponent, 

Rg Reynolds number based on momentum thickness 

T temperature 

u velocity in streamwise direction 

y normal coordinate 

ratio of specific heats ( = 1.4 for air) 

^ boundary layer total thickness 

6* boundary layer displacement thickness 

9 boundary layer momentum thickness 
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p 


density 

Prandtl number {0.72 for air) 


cr 

Subscripts: 
aw adiabatic wall 

c calculated 

e boundary layer edge 

i incompressible 

m measured 

t tabulated 

w wal 1 

DEVELOPMENT OF TURBULENT SHAPE FACTOR EQUATION 


In incompressible flow the ratio of displacement thickness 6*^ to 
momentum thickness 9^ is strictly a function of the shape of the velocity 
profile and is, indeed, called the shape factor. Hence it is possible 
to relate 6*^ and 0^ to N, as described below. 

The incompressible-flow forms of the displacement and momentum 
thicknesses are 



' /'-ST (' - 

® -I 

Assuming a power law velocity profile in the form of equation (1) permits 
the closed-form integration of equations (2). Thus 


0 
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6 *. = 

^ N + 1 


^ (N + 1) (N + 2) ^ ' 

lividing equation (3) by equation (4) yields the incompressible shape 


In compressible flow the situation is not so straightforward since 
the displacement and momentum thicknesses 




Pe 


are no longer strictly functions of the velocity profile shape. The 
compressibility effects contained in the density profile result in 
integrals which do not have general closed-form solutions. The functional 
relationship of equations (6) may be expressed as H = H (N, M0, ) 

For the zero heat transfer case (T^ = T^^^) H reduces to 
(N, Me). Reference 1 has numerically integrated equations (6) for this 
zero heat transfer case, and has tabulated as a function of N and Me 
This function would be much more useful if it could be put in analytic 
rather than tabular form. Hence, an attempt was made in this study to 
analytically define the functional relationship Haw " ^aw ^e)» 
described below. 

Figure 2 was prepared from the tables of reference 1 to show the 
variation of Haw with N and Me • If Hj^i, is viewed as being composed of 
Hi (N) plus a compressibility correction term aH (N, Me), we can write 

Haw = Hi (N) + AH (N, Mg) (7) 
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a 


The variation of Ha^ with Mg appears exponential. Hence if we assume that 
H = Haw - Hi " , the exponent x can be obtained from the slope 

of the logarithmic plot of Haw ^e» tHat is 


aloqo(AH) 

^ ' alogelMe) 

Figure 3 shows that the exponential variation fits the data very well, 
and that the value of the exponent is about 1.989 . 

Hence ^H= f(N) Mg^‘^®^ (8) 


or 


M, 


AH 

1.989 


f(N) 


Figure 4 shows how the parameter 


'AH 


Me 


g - g - j varies with N. The 


function f(N) was estimated from this figure to be 

0.2719 


f(N) = 0.4099 + 


N 


(9) 


Hence the desired relationship between H, N,and Mg for turbulent 
adaiabatic flow is obtained by inserting equations (5), {8)^and (9) into 
equation (7) . Thus 


Haw = 1.0 + 


¥ *(»• 


4099 + 


0.2719 


) Mp^ • 


989 (10) 


For given values of N and Mg, comparisons were made between the 
values of Haw calculated by equation (10) and Tucker's tabulated values. 

For the range 5^N<11 and OsMesS , which is the range on which equation 

< 

(10) was derived, the disagreement between calculated and tabulated Hgyj 
values is less than 0.3 percent. Even up to Mg = 10, which is the limit 
of the tabulated data, the disagreement is less than 0.5 percent. Thus 
equation (10) provides a satisfactory approximation for the turbulent ^ 
adiabatic shape factor as expressed in Tucker's table. 


Persh and Lee ( reference 2) have numerically integrated equations 
(6) for the heat transfer case, and have tabulated H as a function of N, 
M and fw ~ 1^aw ^ In order to get an analytical expression for this 

e > T 

'e 
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general cascj H was viewed in this paper as being composed of the zero heat 
transfer value plus a correction term due to heat transfer SR . Thus 
we may write 


H (n. Me, (N, Mg, 0) + SR Mg, ) 

Figure 5 was prepared from the tables of reference 2 to show the 

variation of AH (=H - Haw) with the heat transfer parameter Iwjl_2.^w 

Te 

for Mach numbers of 0 and 5, and N values of 5 and 11 - This variation 
is approximately linear. Hence we may write 

AH = g(N, Me) (12) 

'e 

Morenver, Mach number has a very small effect on the function g, as shown 
in Table I where the g values were measured from figure 5 . 

TABLE I - MEASURED VALUES OF g 


Fn 

Me 

g 

5 

0 

1.2861 

1.27 

5 

5 

1.256) 


11 

0 

i.ust 

1.13 

11 

5 

1 . 124 ) 



If the dependency of g on Mg is eliminated, eq. 
AH = g(N) - '''aw 


(12) can be written 
(13) 


Assuming a linear variation of g with N, the following equation can be 
derived from Table I . 

g(N) = 1.39 - 0.024N (14) 

Hence equation (13) becomes 

= (1.39 - 0.024N) ~ 05) 

Te 

Comparisons were made between the values of AH calculated from equation 
(15) and those taken from the tables of reference 2 . The disagreement 
between the calculated and tabulated values over the range of variables 
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S 10) 


covered by the tables (0^Me^2O, S^N^Il ,, and - 10 ^ ~ 
was less than ± 2 percent . ® 


Hence the general expression for the turbulent shape factor as a 
function of power-law exponent, Mach number, and heat transfer can be 
obtained by inserting equations (10) and (15) into equation (11) . 


H = 1.0 + 


2.0 

N 


+ ^0.4099 + 


Me 


1.989 


+ (1 .39 - 0.024N) 



(16) 


COMPARISON OF LAMINAR AND TURBULENT 
SHAPE FACTOR EQUATIONS 


An approximate relation for the effect of compressibility on 
laminar shape factor has been derived by Monaghan in reference 3. The 
relation is 


2jrA-jrD-4B-4 
^ 4 -7T 

where A = Li 


Te 



(17) 


Hence equation (17) becomes 




P / \ - Taw 

Mg2 . 4 + [2n - 4 Tp 

4 - 7T ' 


The adiabatic wall temperature ratio in laminar flow is 
Taw s 1 + Me^ 

i 2 


(18) 
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Inserting equation (19) into equation (18) results in the following 
expression for Y= 1.4 and a= 0.72 . 


H = 2.660 + 0.715 + 3.143 



( 20 ) 


Figure 6 shows comparisons between the laminar and turbulent shape 
factors, equation (20) and (16) , respectively. The general trends are 
similar; however, both Mach number and wall temperature have somewhat 
larger effects on the laminar shape factor. For experimental data in 
which the integral parameters 6 * and 9 are known, figure 6 can be used as a 
convenient criterion to determine if a profile is laminar or turbulent. 


CALCULATION OF POWER-LAW EXPONENT FROM SHAPE FACTOR 


Note that the adiabatic turbulent shape factor equation developed 
in this paper - equation ( 10 ) - may be easily solved for the power law 
exponent N. 


N 


^ 2.0 + 0.2719 Mp 


1.989 


1.989 


( 21 ) 


Haw - 1*0 - 0.4099 Mg' 

Thus it should be possible for a given freestream Mach number to calculate 
the N values of turbulent profiles directly from experimental integral 
parameters. 

If the incompressible turbulent shape factor equation (equation (5)) 
is solved for N, the result is 


N = 


H-i 


1 


( 22 ) 


Note that for the case of Mg = 0, equation (21) reduces to equation (22) . 

The calculation of N by this technique can be accomplished by one 
of the follov.'ing methods: (1) For profiles in which « T^^ ^ \\ can be 

calculated directly from Haw ^nd Me by equation (21); or (2) For profiles 
obtained under heat transfer conditions, N can be interpolated from equation 
(16) from known values of H, Me , and Tw - Taw ; or (3) For any profile 
in which the incompressible, or kinematic,^^onns of displacement and momentum 
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thicknesses are known, N can be calculated from H-j by equation (22) . Note 
that this third method is not restricted to Mg = 0 or profiles even 

though Hi is used in the calculations . 

Since integration of experimental data to obtain <5* and 0 is a normal 
data reduction procedure, it would be a simple matter to integrate the same 
profiles to obtain 6*^- and 9^ . Hence by using the appropriate method the 
experimenter can calculate the N values of his turbulent profiles without 
the need of data plotting, curve fitting and slope measuring. 

EXAMPLES OF USE 


Let us now examine a few adiabatic velocity profiles to see how the N 
values calculated by this technique (Nq and N^, i) compare with measured 
values (N„,) . Reference 4 contains 12 turbulent velocity profiles at Mach 
numbers of 1.975, 2.320, and 4.630 . Table II lists the test conditions 
of these profiles, the integral parameters including Nc and Ng^ .j , and 
measured in the conventional manner. 
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TABLE II - REFERENCE 4 DATA 


PROFILE 

Me 


Haw 

Hi 

Nc 

Hc,i 


1 

1.975 

3.62,c104 

3.002 

1.275 

7.35 

7.27 

7.41 

2 

2.320 

1.48x10'^ 

3.594 

1.238 

8.46 

8.43 

8.50 

3 



3.575 

1.228 

8.88 

8.78 

8.63 

4 



3.554 

1.217 

miiQiii 

9.21 

9.06 

5 



3.534 

1.207 


9.63 

9.45 

6 

2.320 

8.04x10^^ 

3.532 

1.204 

9.98 

9.33 

B 

7 

4.630 

1.07/10^ 

10.508 

1.237 

8.91 

8.44 

9.00 

8 


1.69^10^ 


1.233 

9.10 

8.58 

9.20 

9 

4.630 


10.418 

1.214 

9.94 

9.37 

10.20 

10 

. 4.630 

4.44a10^ 

10.369 

1.200 

10.60 

9.98 

10.83 

11 

4.630 


10.343 

1.194 


omul 

11.11 

12 

4.630 

EQSulfli 

10.322 

1.188 



11.31 
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A comparison of the experimental H values of this table with the curves 
of figure 6(a) reveals that these profiles are indeed turbulent. The ratios 
of calculated-to-measured N values are plotted in figure 7 to show that both 
Nc and Nc, i are generally in good agreement with the measured values of N. 
Reference 5 has examined a large number of profiles and has estimated that N 
can be measured within an accuracy of about +10 percent. For all the data 
presented in figure 7> the agreement between calculated and measured N values 
is within this accuracy. Also, note from Table II that at each Mach number 
both the calculated and measured N values show the expected trend of 
increasing N with increasing Reynolds number. 

One of the profiles from reference 4 is examined in more detail in 
figure 8 . Straight lines were fitted to the outer part of the profile only 
(solid line) and the inner part only (dotted line) to illustrate the range 
of curve fits which could be drawn through this profile. The N values of 
these lines are 6.69 and 8.68 , respectively; which represent a ±13 percent 
range. Since it is unlikely that these extreme fits to the data would be 
performed, especially the inner profile fit, the range of N curves which 
would normally be fitted to this profile is probably close to the ±10 percent 
mentioned earlier . A more reasonable fit to the complete profile is 
represented by the dashed line, whose N value is 7.41 , which is very close 
to the calculated values for this profile (Nc = 7.35 and Nc, i = 7.27). 

Since this technique calculates N from either H, Haw« ^^i » which are 
obtained by integration of the entire profile, the N values thus obtained 
represent a fit to the entire profile. Good results would not be obtained, 
therefore, for turbulent profiles containing a large laminar sublayer. 
Measurements in a laminar sublayer tend to reduce the calculated N values, 
with the extreme example being a completely laminar profile. For example, 
the integration of a Blasius profile yields Nc, i =1.26 . Since N ranges 
from about 5 to 11 for turbulent profiles, this technique could also be 
used as convenient criterion for determining whether a profile is laminar or 
turbulent. 

One cautionary note about this technique needs mentioning. It can be 
determined from equation (21) that N is rather sensitive to errors in H. 

For example, integration errors which result in a one percent error in H 
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produce about a 10 percent error in N . Therefore in order to use this 
technique, accurate integration of experimental data is required. 


CONCLUDING REMARKS 


In summary, analytical expressions for the effects of compressibility 
and heat transfer on laminar and turbulent shape factors H have been 
developed. Solving the turbulent equation for the power law velocity profile 
exponent N has resulted in a simple technique by which the N values of 
experimental turbulent profiles can be caulculated directly from the integral 
parameters. Thus the data plotting, curve fitting, and slope measuring, 
which is the normal technique of obtaining experimental N values, is 
eliminated. The N values obtained by this method should be within the 
accuracy with which they could be measured. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Virginia, August 20, 1974 
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Figure 2.- Mach Number Effects on Adiabatic Turbulent Shape Factor 



d log. (AH) 
s 

d log. (Me) 


1.989 


toge (Me) 


igarithmic Variation of aH with Mach Number 


.50 



Figure 4.- Variation ot Compressibility Parameter with Power-Law Exponent 

















